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Abstract 

The complexity of a neuronal cell shape is known to be related to its 
function. Specifically, among other indicators, a decreased complexity in 
the dendritic trees of cortical pyramidal neurons has been associated with 
mental retardation. In this paper we develop a procedure to address the 
characterization of morphological changes induced in cultured neurons by 
over-expressing a gene involved in mental retardation. Measures associ- 
ated with the multiscale connectivity, an additive image functional, are 
found to give a reasonable separation criterion between two categories of 
cells. One category consists of a control group and two transfected groups 
of neurons, and the other, a class of cat ganglionary cells. The reported 
framework also identified a trend towards lower complexity in one of the 
transfected groups. Such results establish the suggested measures as an 
effective descriptors of cell shape. 



1 Introduction 

Classical investigations in neuroscience have added a large amount of descrip- 
tive data to the knowledge of the mammalian neural system and ultimately of 
the functioning of the brain. In recent years the use of computer and algorith- 
mical methods has been incorporated into the biological sciences with the aim 
of achieving quantitative analysis and modelling. Being no exception to this 
important tendency, mathematic-computational investigations in neuroscience 
are still in an incipient stage, with many challenging questions and potential 
for important results. While the effort to model quantitatively the emergent 
features of a neural system or element (e.g. a single cell) involves its own tech- 
nical difficulties, neuroinformatics has already led to important results regarding 
the interplay between the morphology of dendrites and their firing patterns pQ. 
One difficult aspect implied by the use of computers in neuroanatomy is the 
relatively limited availability of accurate spatial data of real dendritic arboriza- 
tions for diverse morphological classes. The alternative approach of generating 
morphologically realistic neurons has been successfully carried out by using a 
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complementary approach based either on recursive |2] or topological E] algo- 
rithmic implementation. 

The anatomic details of the neuronal cells have been shown to represent a rel- 
evant feature, for phylogenic and ontogenic studies, in the diagnosis of diseases 
and in investigating the interplay between form and function 01210]. With 
the availability on the Internet of combined databases of virtual and real (histo- 
chemically marked) neurons, the analysis or morphological characterization of 
newly acquired neuronal information can now be the subject of systematic pro- 
cessing by computer methods, provided suitable shape descriptors are identified 
and applied. 

There is a current interest in many branches of Physical and, more recently, 
Biological sciences to study in a quantitative way the geometrical outline of 
structures, both static and in development, and from that knowledge to probe 
further into the functioning and justification of those forms. Examples of results 
of such a quest in biology ranges from the elucidation of the Fibonacci spirals 
emerging in sunflowers [7] [5j to the fractal characterization of the nature of 
many botanical plants [jjj. For over a century the search aimed at establishing 
a relationship between the form and function of neuronal cells has challenged 
and stimulated researchers, starting from the pioneering efforts of Ramon y 
Cajal ^UJ. While many diseases have been diagnosed based on the visual per- 
ception of morphological alterations in general tissue morphology, automated 
approaches remain promising subjects of research, with many potentially im- 
portant applications. 

Although it is now generally accepted that dendritic morphology plays a cru- 
cial role in the functioning of the neural cell and ultimately in the behaviour of 
the neural system ^2 H2I El E] > there is no deterministic way to choose the 
best quantitative descriptor of the geometry or topology of single neurons. In or- 
der to be particularly useful, such descriptors should correlate in some way with 
the respective performed functions and behaviour, allowing abnormalities and 
specific divergences from a healthy state or development to be clearly identified. 
An additional benefit of such an investigation of the physical basis of neuronal 
shape and growth is the possibility for investigating, through modelling, the 
relationship between neuronal shape and function. 

The study of biological forms requires the selection of shape descriptors or 
shape functional that fulfil a set of requests. First, measures are expected to 
be objective and fast (in order to treat a representative number of cases), while 
retaining comprehensiveness and being potentially discriminative. In addition, 
such measures should lead to meaningful biological interpretations. An impor- 
tant theoretical issue is the relative degeneracy of a set of measures, a concept 
that reflects the fact that many different forms or shapes produces the same 
or approximated measures. It has been experimentally shown |51 1151 Ho] that 
multiscale analysis tends to augment the resolution power of geometric descrip- 
tors and, despite the increased computational demands, this procedure provides 
useful additional information and characterization of the considered objects. 
The neuronal cells have, in general, a complex spatial tree-like structure with 
many dendritic bifurcations leading to efficient spacial coverage, amplifying the 
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influence area and optimally connecting the neuron with a neighbouring cell's 
dendrites or axons. This structure has been shown to exhibit different levels of 
partial fractality at different spatial scales El El EH ■ 

A recently reported El E] procedure for calculating additive shape func- 
tionals, known as Minkowski functionals gathered in a framework called Integral- 
Geometry Morphological Image Analysis, MIA, has been successfully applied 
to many areas of research including Statistical Physics, Cosmology and Mate- 
rial Sciences [21 El El 1201 E] • A first attempt to bring those morphological 
concepts into biological sciences, targeting the efficient discrimination of two 
main types of the domestic cat's ganglionary neuronal cells, has been reported 
recently |24| . The use of additive functionals is implemented in a pixelwise ap- 
proach leading to the possibility of straightforwardly obtaining the above men- 
tioned multiscale fractal dimension for complementary shape description. Every 
additive continuous and motion invariant functionals in the 2D Euclidean plane 
can be expressed as a combination of three Minkowski functionals which are pro- 
portional to a known geometric quantity, namely the metric area and perimeter 
and the topological Euler connectivity number. This completeness extends to 
higher dimensions and the geometrical functionals multiplicity is always one 
plus the dimension of the lattice generalized voxels (e.g 3 in 2D plane, 4 in 3D 
space, etc.). While the theory of Integral geometry provides a sophisticated 
set of results and formulae, the practical implementation of the procedures in 
Image analysis is relatively simple and efficient El EH ■ The basic ideas of this 
methodology is described in section El and a simple example is worked out to 
outline the procedure. 

In this paper we investigate the importance of the spatial distribution of 
branching points in neuronal arborizations to provide a discriminative and in- 
formative characterization of the neuronal morphology, with special attention 
given to additive shape functionals. This novel procedure involves mapping the 
neuronal image onto a set of points representing its bifurcation pattern, which 
are subsequently dilated in order to produce a multiscale representation E], as 
the Minkowski functionals of such a set of points is recorded. The obtained re- 
sults substantiate the potential of the above procedure regarding its application 
to a database of 3 categories of rat neuronal cells and one class of cat ganglion 
cells. The results are in accordance with the biological importance of the spa- 
tial distribution of branches ^1 ESI EH El El \U\ and provide a discriminative 
measure for the morphological characterization of the neurons, indicating subtle 
morphological differences as a consequence of the considered gene transfectation 
in the rat categories. This suggests that those treatments may not significantly 
alter the topological aspect of the branching cell form as far as the considered 
measurements are concerned. 

2 Additive Functionals 

Integral geometry algorithms have been successfully used to characterize mor- 
phologically complex patterns where the precise process of formation is not 
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completely known and is a subject of modelling, see for example [2BI - The cen- 
tral procedure is the calculation of intrinsic volumes or Minkowski functionals 
(or yet querrmassintegrals) , which are a generalization of the usual determina- 
tion of volume. They can be defined as related to both differential and integral 
geometry setups. 



Figure 1: The disjoint decomposition of a pattern K (left) into a sum V (right) 
of disjoint interior elements. 

For the sake of a better explanation, we start with a description of the 
practical aspect of the adopted procedure. The Minkowski functionals of a 
body K in the plane are proportional to its area A(K), perimeter U(K) and 
the connectivity or Euler number x(K). The usual definition of the connectivity 
from algebraic topology in two dimensions is the difference between the number 
of connected n c components and the number of holes rih, 

x{K)=n c -n h . (1) 

So, if K happens to be an image of the number 8 as in Figure ^ its con- 
nectivity number would be x{K) = 1 — 2 = — 1. One instance where these 
functionals appear naturally is while attempting to describe the change in area 
as the domain K, now assumed to be convex, undergoes a dilation through a 
parallel set process using a ball B r of radius r 

A(K®B r ) = A{K) + U(K)r + Trr 2 . (2) 

Generalizing to higher dimensions, the change in hyper volume is given by 
the Steiner formula [2U 

v d (K®B r ) = J2( d )wi d) (KK, (3) 
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where the coefficients Wv are referred to as Minkowski functionals. For 
instance in the plane (d = 2), 

WP(K) = A{K)M 2 \K) = W., W ^(K) = X (K)n. (4) 

Despite the wealth of results and continuum formulae for obtaining these 
functionals, it is useful to recur to the discrete nature of the binary images we 
wish to analyse by looking at the distribution of pixels in the planar lattice. 
On that level, by exploring the additivity of the Minkowski functionals, their 
determination reduces to counting the multiplicity of basic building blocks that 
compose the object in a disjointed fashion. Figure ^ illustrates this process. A 
pixel can be decomposed as a disjointed set of 4 vertices, 4 edges and one square. 
The same process can be applied to any object in a lattice. The fundamental 
information needed here is a relationship for the functionals of an open interior 
of a n-dimensional body K which is embedded into a d-dimensional space 

W^(K) = (-l) d+n+v wW(K),v = 0, . . . ,d. (5) 

As there is no overlap between these building blocks and by the property 
of additivity of these functionals, for a body P composed of disjointed convex 
interior pieces N m , we may write 

WW(V) = J2w!) d) (N m )n m (V), v = Q,...,d. (6) 

m 

Where n m (V) stands for the number of building elements of each type m 
occurring in the pattern V . For a two-dimensional space, which is our interest 
regarding the considered neuronal images, we display in Table ^ the value of 
Minkowski functionals for the building elements in a square lattice of pixels and 
their direct relation to familiar geometric quantities on the plane. Using the 
information (with a = 1) presented in Tableland equation JJjJ we have 

A(V) = «2, U{V) = -4n 2 + 2m, X (V) = n 2 - n x + n Q . (7) 

Going back to Figure ^ we find, for this specific example, as n 2 = 16, 
ni = 47 and no = 30, that A = 16, U = 30 and X = — 1. So the procedure to 
calculate Minkowski functionals of a pattern K has been reduced to the proper 
counting of the number of elementary bodies of each type that compose a pixel 
(squares, edges and vertices) involved in the make up of V . 

3 Branching point patterns 

The existence of profuse branching structures in natural shapes provides ev- 
idence of the effectiveness of such shapes in providing an interface with its 
environment as well as for enabling connections. As pointed out in |12j many 
conditions may interfere with the morphology of the neuronal tree structure, 
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Table 1: Minkowski functional of open bodies N m which compose a pixel K: 
P (vertex ), L (open edge) and Q (open square). 
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for instance learning, 'enriched' environment, hormonal fluctuations and levels 
of bioelectric activity. 

The relationship between dendritic morphology and cell functioning, alone 
or connected in networks, has drawn the attention of scientists leading to a 
search for a set of measures that would as completely as possible describe the 
neuronal shape [221 • For example, the dependence of the dendritic diameter on 
the distance from soma, the relationship between the dendritic diameter before 
a bifurcation point and the diameter of the two daughters stemming from this 
point, or the ratio of diameters of daughter dendrites. These are standard 
examples of a global feature (the first) and two local properties (the following 
two) one can specify/calculate in describing/analysing a neuronal shape, real 
or virtual, see |2J. There is a vast recent literature describing shape analysis in 
general and applied specifically to neuronal shape, see for example 1301 EH 
for a general dendritic description, and ^H] for multiscale fractality, ^5] f° r 
wavelets and [S] for bending energy applied to ncuromorphometry. 

Since the pioneering work of Sholl > some authors have mentioned specific 
situations where the well-known Sholl analysis could lead to degenerate results, 
assigning to visually different neurons the same (in statistical terms) descriptor 
values, see [SU] for example. The method makes use of a reference point at the 
soma and draws concentric circles through the dendritic field, while counting 
the number of intersections within each circle. While this process is easy to 
implement manually, which may account for its popularity, it clearly fails for 
asymmetric forms. While alternative approaches |31| report success in improv- 
ing this setup, we take a different route, abandoning the soma as a reference 
point. 

We start our analysis by taking all bifurcation points and constructing an- 
other image keeping the metric relation between these points. Although there 
are automatic ways of extracting salient points |221ES1I^; we implemented this 
task in a semi-automated way, which is followed by a sequence of exact dilations 
(permitted parallel set dilation ,[35]) of the selected points. Although primar- 
ily interested in the connectivity, all Minkowski functionals are calculated by a 
routine that counts efficiently the multiplicity of the building elements for each 
neuron image at every radius of dilation as described in section [3 equation [7| 
The idea is to use the connectivity of the set of points to capture, at each scale, 
the spatial relationship of the neuron branching structure. As this measure is 
invariant to rotation and translation, as well as to scaling, eventual variations 
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of size due to the relative stage of development can be disregarded. 

We investigate the morphology of two different categories of cells. One in- 
volves cat ganglion alpha cells, which constitutes a class by itself but rather 
diversified in form, see [351 124) . The other includes foetal rat cerebral cortex 
neurons cultured in a 2D tissue system This category is subdivided into 
three differently treated classes consisting of a control group, a group of neurons 
transfected with a gene OPHN1 which encodes for oligo-phrenin, and a positive 
control group transfected with pl90 RhoGAP. 

Figure |21 shows examples of two treated rat neuronal cells, illustrating the 
branching selection process adopted in this work. Note that the neuron images 
have been rescaled in this picture, but the final dilation radius in all processed 
samples are the same and equal to 30 pixels. 

It is important to comment on the complementary nature of the framework 
proposed in this article with respect to more traditional approaches reported 
in the literature. When compared to the previous Sholl methodologies, the 
measurements in our approach adopt multiscale concepts so as to retain as 
much as possible the metric and global structure, while avoiding the somewhat 
arbitrary use of the soma as a general reference. In order to express in a more 
objective and effective way the morphological features of special interest, namely 
the overall complexity of the neuronal pattern, our procedure was designed from 
the start to concentrate on the bifurcation points, which are used as primary 
data. No attempt is made here to achieve completeness in the shape description, 
or to use the considered measures as features to be incorporated in an algorithm 
to grow artificial neurons. However, the reported measurements can be useful for 
statistical validation of such artificial data, in the sense that both the simulated 
and original neuronal structures should lead to the same probability densities. 

4 Results 

For each cell we calculate multiscale curves of connectivity as shown in Figure 
As discriminating measurements we take the area under the interpolating curve 
(its integral) and the standard deviation of the difference between that curve and 
the original data points. The former can provide information about the overall 
structure, while the latter characterized the finer details. Having performed 
similar analysis for area and perimeter functionals, as well as for the multiscale 
fractality as derived from the area data, we decided to focus on the multiscale 
connectivity of the branching points set as a descriptor for its relationship to 
the spatial distribution of branching points and its above average discriminating 
power. 

Our primary concern here is to investigate the effectiveness of the procedure 
adopted to quantify the subtle morphological aspects that may characterize 
different cell types while also expressing a biologically relevant attribution of 
the cell form, namely the distribution of bifurcation points. In this regard, our 
results, shown by Figure |3J led to a well-defined separation of the alpha cells in 
a feature space defined by the integral of the multi scale connectivity and its 
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Figure 2: Isolating morphological information relating to the neuronal branching 
structure while preserving the metric relationship among the branching points. 
Each bifurcating point undergoes an gxact dilation while its connectivity is 
recorded as function of the scale parameter. Prototypical transfccted rat neu- 
ronal cell, a) and d). Extracted pattern of bifurcation points, b) and e). The end 
of the bifurcation point parallel set dilation procedure, c) and f). The resulting 
discriminating functional for both neuron types, g). 
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Figure 3: A discriminative feature space. Most of the Rho-Gap transfected 
neurons went to the lower left corner of the plot, while the Alpha cat cells are 
spread by themselves towards higher values. There is a fuzzy display of both 
control and OPNH-1 samples. 
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Figure 4: The discriminative power of the above feature space in depth. A 
discriminant analysis reveals a reasonable separation among ALPHA and Rho- 
GAP groups and is not conclusive for the others. 
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Table 2: Canonical Coefficients 
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Table 3: Mahalanobis Distance 
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standard deviation. 

Although the three differently treated rat cells are not clearly distinguish- 
able, one group of them, namely the one containing neurons which show over 
expression of RHO-GAP following transfectation, clearly tends to populate the 
lower left corner of the feature space. Yet, while high dispersion is shown by 
the rats sample cells, including the control group which is as much diversified in 
form as the treated groups, the statistical tendency towards complexity agrees 
with a previous analysis which focused on the multiscale fractal dimension of 
this group of cells [3*7] . 

Figure [""J shows a scatter plot defined by two canonical discriminant func- 
tions (their canonical coefficients are displayed in Table [""J for the above se- 
lected shape descriptors. The separability in groups is optimized (see Table 01 
for a quantitative index, the Mahalanobis distance) in such an analysis, leading 
to a somewhat improved visibility of the above mentioned distinction of the 
Rho-GAP over expressing set of cells from the others. Nonetheless the other 
differently treated group mixes up strongly with the control group, as can be 
verified from Table 21 which presents the plug- in (also called confusion matrix) 
classification results. As can be seen from this statistical analysis the alpha cat 
cell is not trivially distinguishable from the others, as shown by the overlaps 
in both the feature and canonical discriminant scatter plots. While the effec- 
tiveness of use of Minkowski functionals (considering the whole cell) has been 
previously reported |24|. the framework developed here is meant to emphasize 
the distribution of the bifurcation points, in a way similar to that in Sholl's 
analysis, on the cell morphology. 

5 Conclusions 

This work describes how promising results have been obtained regarding a novel 
procedure and measurements for classifying different types of neurons and to 
reveal morphologically relevant attributes of different classes. Although the 
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Tabic 4: Plug-in classification table 
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subtle morphological variations induced by the treatment of our samples are 
masked by the statistical distributions and show only a small tendency towards 
a reduction in complexity, the procedure shows clearly its discriminating power 
when applied to a well characterized class of alpha cat ganglionary cells. As the 
considered measurements were designed to express a particular biological aspect 
of the neuronal geometry, namely the spatial relationship among the branch- 
ing points, we would expect that this specific trait is not significantly affected 
by the gene transfectation process. The procedure reported in this paper can 
be generalised to process 3D Neurons, real or virtual, and can be implemented 
straightforwardly. By extending the present analysis to 3D we suspect that sub- 
tle information relating to neuron complexity could be significantly improved. 
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